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Two leading dark matter candidates from supersymmetry and other theories of physics beyond the 
standard model are WIMPs and weak scale gravitinos. If the lightest stable particle is a gravitino, 
then a WIMP will decay into it with a natural lifetime of order a month ~ M^i/M^^^j^. We show 
that if the bulk of dark matter today came from decays of neutral particles with lifetimes of order 
a year or smaller, then it could lead to a reduction in the amount of small scale substructure, less 
concentrated halos and constant density cores in the smallest mass halos. Such beneficial effects 
may therefore be realized naturally, as discussed by Cembranos, Feng, Rajaraman, and Takayama, 
in the case of supersymmetry. 



I. INTRODUCTION 

We have many independent lines of observational ev- 
idence that point to the existence of dark matter in the 
universe. The dark matter particle, however, has not 
been observed. Most theories of physics beyond the stan- 
dard model that have new physics at the weak scale, nat- 
urally predict the existence of WIMP (Weakly Interact- 
ing Massive Particle) dark matter. The weak interac- 
tion cross-section endows the WIMP with the right relic 
density to make up a substantial fraction of the energy 
density of the universe. Thus most cosmological studies 
of dark matter assume that it is a stable cold particle 
with interaction cross-sections comparable to the weak 
interaction cross-sections. 

Within theories of physics beyond the standard model, 
one also has the gravitino (or other gravitationally cou- 
pled particles). Recent work by Feng, Rajaraman and 
Takayama provided an interesting twist to the old 
dark matter tale. They observed that if the gravitino 
is the lightest supersymmteric particle (LSP), then the 
next lightest particle (NLSP) can decay into the grav- 
itino with a lifetime that is long since the gravitino cou- 
ples only gravitationally. The natural time scale for this 
decay is M^^/M^^^^^ ~ month. This opens up the pos- 
sibility that at least a fraction of the dark matter we 
observe today comes from decays when the universe was 
days to years old. We note that similar models may be 
constructed in the context of cosmologies with extra di- 
mensions 01 ■ Two early universe constraints on these 
models are that the electromagnetic decay products must 
not distort the black body spectrum appreciably and the 
light element abundances from BBN must match with 
observations. Recent work has shown that there are re- 
gions of parameter space where all these conditions may 
be met 

II II H II 0- In the present work we work out 
the late time consequences of these decays. 

Another way to get dark matter to decay around sim- 
ilar epochs is to kinematically suppress the decays such 
that the lifetime is orders of magnitude larger than the 
"natural" time scale. Profumo et al. Q provided the 
first example of such a decay that arises in supersymmet- 
ric theories. The consider the decay of stau into the LSP 



with a mass splitting smaller than the tau lepton mass. 
In this case, the two body decay is forbidden and the stau 
can only decay into 4 daughter particles. The 4-body de- 
cay of a charged slepton is suppressed by the NLSP-LSP 
mass splitting as (M/AM)^. Thus with AM/M - lO'^, 
the stau lifetime is of order a few weeks Q . 

Sigurdson and Kamionkowski and then Profumo et 
al. Q worked out the impact of the charged stau decay on 
the matter and CMB power spectrum. They found that 
the dark matter power spectrum was suppressed due to 
the coupling of the charged NLSP to the photon-baryon 
system. This coupling damps the matter power spectrum 
once a perturbation mode enters the horizon. 

In this paper we point out that the large suppression of 
the matter power spectrum is possible even if the NLSP 
is not charged. Unless the mass splitting between the 
NLSP and LSP is fine-tuned to be small, we point out 
that the decay will impart a large velocity to the LSP. 
Effectively the LSP will act like warm dark matter. This 
leads to a cut-off in the matter power spectrum. We 
calculate shape of the power spectrum given the lifetime 
and the masses of the particles involved. 

We also calculate the effect of the smaller phase space 
density that results from this decay on dark matter halos. 
They could (for small enough primordial phase space den- 
sity) be instrumental in shaping the smallest mass halos, 
lowering the fraction of mass in sub-halos, and reducing 
the concentration of larger mass halos. Stringent con- 
straints on these decays come from the requirement that 
the universe be reionized by a redshift of six or higher. If 
the present hints [T9 | for early reionization are borne out 
by future data, then this could rule out a big chunk of the 
supersymmetric parameter space in which the gravitino 
is the lightest particle. 

Cembranos et al. ,Xl] show that the beneficial effects 
mentioned above may be obtained in a large region of 
the super-symmetric parameter space. The prospects for 
testing this region of parameter space wherein we have 
early decays is promising. These models will be tested by 
observations made with NGST of the high redshift uni- 
verse (e.g., (.12,|). by future surveys of strong lens systems 
(e.g., [13), weak lensing measurements (e.g., |lj|), as 
well as the more traditional CMB and large scale struc- 
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ture observations. 

The dark matter observed today could also be a mix- 
ture of particles produced during the reheating phase and 
those arising from the decays. For the gravitino example, 
the two contributions are expected to be comparable for 
a reheating temperature of 10^" GeV. This kind of mixed 
model would ameliorate the rapid drop in power on small 
scales. 

Dark matter from decays have a non-thermal momen- 
tum distribution. Lin et al. [l5| . Hisano et al. |l6j| . 
Kitano and Low 0| studied different models of non- 
thermal WIMP dark matter and their effect on the small 
scale structure. The generic cosmological effects they 
found are similar to what we discuss here because the 
effects stem from the large velocity dispersion. Decays 
leading to a large velocity dispersion were also consid- 
ered briefly by Hogan and Dalcanton |0| in their work 
on the astrophysical consequences of a small primordial 
phase space density. 



II. DARK MATTER FROM DECAYS 

Consider the decay process, DDM DM -I- L, where 
L denotes "light" particles with mass mi , DM is the dark 
matter today with mass m and DDM is the parent par- 
ticle with mass M. When we derive the perturbations, 
we will be enforce the limit mi << m as is appropri- 
ate for supersymmetric models (e.g., 0). We will also 
assume that DDM and DM are neutral, for example, 
a sneutrino decaying into a neutrino and a gravitino. 
All of the effects we mention here are relevant also for 
charged particle decay. In addition, the coupling of the 
charged parent particle to the photon-baryon system re- 
sults in the damping of perturbations that come into the 
horizon when the charged parent particles dominate the 
non-relativistic matter density |^. We will work out the 
physics and astrophysics of charged particle decay in a 
future publication. 

For a generic n-body decay, neglecting Pauli-blocking 
factors and inverse decays, the change in the phase space 
due to decays is given by 



/ddm(pddm) 
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(1) 



where -Eddm = Pddm + ^^-d /ddm is the phase space 
distribution of the DDM particles. The over-dot denotes 
derivative with respect to the conformal time, drj = dt/a 
and a is the scale factor. Specializing to two-body decays, 
one can show that the DM phase space is populated by 
the DDM decays according to the equation 
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where ^2,1 = iO-^Euuml ± pBMPcMM)/ml,^. pcM is 
the center-of-mass momentum and ttiq = Af^ -I- m^y^ — 
ttil . An analogous equation holds for the other daughter 



particle. Note that for a two-body decay we have life- 
time T — SttA/P /pcmIMI"^ where |A^p is the quantum 
mechanical amplitude for the decay process. Also note 
that since the decay happens when the DDM particle 
is cold, PCM is an accurate measure of the momentum 
imparted to the daughter particles. 

The equations for the change in the phase space dis- 
tribution due to the decays, Eq. ^ and Eq. [21 can be 
integrated to yield the following equations for the change 
in the density of DM and DDM. 



PDDM + 3-/7DDM 



Pdm + 3-(pdm + Pdm) 
a 



-PDDM 



(3) 
(4) 



where Pdm is the pressure of DM and it is comparable 
to Pdm at early times when the bulk of decay occurs. 
An equation for "L" particles may be written down by 
inspection of the DM equation. 

We will also need to calculate the phase space distri- 
bution of DM particles. One can show that this is (in the 
limit of completely non-relativistic decay) given by 



mq^T 



exp(-^)e(apcM - q) 

T 



/dm (9, a) 
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(5) 

where q is the comoving momentum of the DM particle 
and tq = t(a = q/pcAi)- We have assumed that the 
decays are happening during radiation domination and 
hence tq cc q^ and /dm(<7) q^^ cxp{—q'^ / [2p'^]^jTd]) . 



III. PERTURBATIONS 

We now write down the equations for the perturbations 
in DDM, DM and L particles. The first result is that the 
perturbation in density relative to the mean for the DDM 
particles is unchanged. This is derived by using the fact 
that the DDM particles are non-relativistic. The physical 
content of this statement is that since the DDM particles 
do not have a large peculiar velocity during decay, the 
only effect of the decay is to remove the same fraction of 
the DDM particles from every region of space. 

To calculate the perturbations explicitly, we write it as 
a sum of two terms: 

<5/DM,£(fc, q, a) = (5/^M/(^' 9' a) + '^/dm/(^' 9' a) • (6) 

The first term Sfj^^ ^(fc, q, a) is chosen to satisfy the col- 
lisionless Boltzmann equation for a massive particle and 
may be calculated in a manner analogous to the pertur- 
bations in the massive neutrino phase space 21] . The 
second term is harder to calculate numerically. One can 
make approximations similar to the one made in the pre- 
vious section and obtain: 



fuuiq, a)h{k, aq)ji{kujq{a, Qq)) , (7) 
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where Og = q/pcM-, (^q{x,y) = da{q/e{q,a))/a and 

e{q, a) = m?a^ + q^ is the comoving energy. Note 
that ujq{x,y) is the comoving distance traversed by a 
particle with comoving momentum q between x and y. 
h{k,a) is the perturbation to the trace of the the syn- 
chronous gauge spatial metric in fourier space written as 
(5y (5(k)/3 + 5,jh{k, a) + {kjij - 5pp)a{k, a) [p. 

We have modified CMBfast [221 to incorporate the 
above physics. The physics of how the decay affects the 
dark matter perturbations is simple. The decay of DDM 
to DM provides a large momentum to the DM particle. 
This implies that the mean free path of DM is larger than 
in the standard case. Perturbations on scales smaller 
than the "mean free path" of the DM particle cannot 
survive. 

How is this manifested in the Boltzmann equation? 
The solution of the linearized coUisionless Boltzmann 
equation (valid when the decay term can be neglected) 
can be written as 

fO-o poo ir 

Sp{ao) = {27t)-\^' / A / dieq^^ [ 

Jo Jo dq 

- jjoiku;,{a,0)) + ^a'j2{kLUqia,0))] , (8) 

where h' and a' denote derivatives with respect to scale 
factor a. The free-streaming effect comes in through the 
effect of the function ujq. We note that the jo and j2 terms 
result from the plane wave expansion of exp(ifctjq(a, 0)). 
In the limit of m — > 0, we have a plane wave (which 
is how the photons free-stream after decoupling). In 
the other limit of m — » cxd, jo ^ 1 and j2 0, and 
we obtain the evolution of the CDM density perturba- 
tions. Similar expressions were first derived by Bond 
and Efstathiou and Brandenberger, Kaiser and 

Turok (24J. On small scales, one may make further ap- 
proximations and write Sp^ao) = {27r)^^ f^" th{h' + 
a')/3 / dqq^efod\n{wq)/d\n{q)cos{kujg). Note that when 
ku!q(a, q/pcAi) gets large, we get no contributions to the 
integral due to cancellations from rapid oscillations. The 
power spectrum is thus suppressed. 

A point of interest in this discussion is that 
the quantity that determines the damping of the 
power spectrum on small scales is not the canon- 
ically defined mean free path f^'"' dtv{t) / a{t). 
Damping occurs when kujq{a,q/pcM) > 1; very 
roughly, this condition picks out scales smaller than 
0.005(Q/[Mo/pc3/( km s-^])'^/^ Mpc. 

What are the natural variables to describe this class 
of "warm dark matter" models? We may take a^cc and 
m/pcM to be the two variables, adcc is the scale factor 
of the universe when the age is equal to the lifetime r. 
The cosmological consequences depend only on these two 
variables. 

If we specialize to the case of decays in the radiation 
dominated era, then we will find that there is just one 
variable that adequately describes the cosmological ef- 
fects. For those decays that release EM energy, decays 



must occur in the radiation dominated era in order to not 
distort the CMB black body spectrum beyond the 10"'' 
level. 

A physically motivated variable is phase space density 
of DM particles in the early universe. We will adopt the 
definition of Hogan and Dalcanton [3 who defined Q = 
For times much larger than the decay lifetime, 
we have the exact relation 

Q = 10-24/" y McJpc' 

VPCA/adcc/ ( km s i)3 

We will see that as a first approximation, the power spec- 
trum only depends on Q. For reference, we note that in 
the case of the slepton decaying to the gravitino, Q ~ 
2.1 X 10-^{2pcM /Mf{M/ TeY)^-^ Mq/ pc^ / { km s-^f. 

The net result on the power spectrum can be writ- 
ten in a manner analogous to the massive neutrino case. 
The power-spectrum is suppressed on scales smaller than 
Ac- Here we provide a simple fitting formula that en- 
capsulates the basic features of the power spectrum at 
the 25 % level for lifetimes less than about a year and 
k < 25h/ Mpc. 

Pumik) = PcDM Qcxp(-(fcA,)V2) 

+ l[l + (kX,f]-'y , (10) 

where Ac = 0.0198/Q°-2'^^ Mpc. At this level of accuracy, 
we thus have a one parameter family of models. 

The above power spectrum should be compared to the 
power spectrum in Warm Dark Matter (WDM) mod- 
els. Bode, Ostriker and Turok quote Pvvdm = 
-Fcdm(1 + (afc)^'')"^"/'^ based their analytic work and 
the exact results of Ma |2^, where = 1.2 and 
a = 0.048/1-1 Mpc(OwDM/iV0.169)°-i5(keV/m)i i5 for 

a fcrmionic WDM with two internal degrees of freedom. 
Note that the asymptotic power-law in the two models 
are different; the decay model has comparatively more 
power due to its non-thermal distribution. Looking at 
the variance <j(M), we find that the modified DM power 
spectrum with Ac = 0.1 Mpc/h is a good fit to the vari- 
ance function for a 1 keV WDM model. Bode, Ostriker 
and Turok [2^ find that the large scale structure con- 
straints are met by a WDM model with 1 keV particle. 
Thus, Ac < 0.1 Mpc/h will also reproduce the large scale 
structure we observe. For reference we observe that a keV 
WDM implies Q = 5 x 10'^ Mq/ -pc^ /{ km s''^)^. From 
the approximate formula following Equation 1101 we get 
a similar Q value for Ac — 0.1 Mpc/h and h = 0.65. 

IV. PHASE-SPACE DENSITY CONSTRAINTS 

In warm dark matter models like those due to massive 
neutrinos, an important constraint for late time cosmol- 
ogy results from the finite phase space density [Tsl |2^ . 
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For a warm dark matter fermionic particle which decou- 
ples when relativistic, the phase space density frozen in is 
a Fermi-Dirac distribution. The maximum allowed phase 
space density for a Fermi-Dirac distribution is given by 
h~^/2. The resulting gravitational distortions of the 
phase-space sheets can never exceed this phase space 
limit for coUisionless particles. This is a direct conse- 
quence of Louiville's theorem as applied to collision-less 
systems (also called the Vlasov equation). 

For the decaying dark matter particle at hand, the 
maximum phase space density, though finite, can be very 
large because of the 1 /q term. We note that q tends to 
a constant when the parent particle is relativistic. 

The maximum phase space density argument is not the 
strongest statement that can be made regarding the evo- 
lution of collision-less systems. Consider a phase space 
distribution that is Fermi-Dirac plus a delta-function at 
some small q value. The maximum phase space density 
argument would have nothing to say in this case. How- 
ever, if the total number of particles within the delta- 
function spike is small, then we don't expect them to 
affect the evolution of the collision-less system. 

Lynden-Bell showed that collision-less systems 
have an infinity of conserved quantities. He labeled these 
M(> /), the mass of particles (or phase space cells) with 
densities greater than a value /. There is however a 
problem in using M(> /) to make statements about the 
evolution of a collision-less system. Lets consider the 
case of a uniform density of particles collapsing to form 
a galaxy. The observational information we have about 
the galaxy is never about the fine-grained distribution of 
the system. We only measure averaged quantities - the 
coarse-grained distribution. M{> /) is not conserved 
for coarse-grained distributions; it could increase or de- 
crease. Tremaine et al. proved a theorem using the 
Boltzmann H-functionals that coarse-graining decreases 
these H-functionals (concave functions of the distribution 
function, like the entropy). What about comparing two 
coarse-grained systems? 

Dehnen (see also the early work by Mathur [sj) 
recently proved that there exists a function - the ex- 
cess mass function - that always decreases as a result 
of coarse-graining. The more coarse-grained, the smaller 
this excess mass function. The excess mass function is 
very close in spirit to Lynden-BeU's M{> /). It is defined 
as: 

D{f) = J d'xd\{F{^, q) - /)e(i^(x, q) - /) (11) 

Let us look at the two important features of this func- 
tion in a little more detail. First, for coUisionless systems 
D(f) is an invariant. In a coUisionless system, the phase 
space density is conserved along a world-line (solution 
of the equations of motion). Thus the total "mass" in 
elements with F > f must be the same (even though 
the global distribution of this mass could have changed) . 
Second, for coUisionless systems, £*(/) decreases as a re- 
sult of coarse-graining. In order to see why this is so. 



let us coarse-grain by dividing up our phase space into 
"macro" phase space cells . Each macro-cell of phase 
space is made up of many (at least more than one) micro- 
cells. Suppose that the macro-cell created in this way has 
F > f. Now those micro-cells in this coarse- graining that 
had densities larger than / don't change D{f). However, 
the micro-cells with densities smaller than / that get in- 
cluded in D{f) (post coarse-graining) do change D{f). In 
fact, it can only decrease D{f) because the F — f is neg- 
ative for them. We can similarly argue that a macro-cell 
created with F < f also decreases D{f). Thus coarse- 
graining can only decrease D{f). 

This provides a natural way to constrain the inner 
cores of dark matter halos. D{f) for the halo must be 
smaller than the D{f) calculated for the primordial dis- 
tribution function for all values of /. We note that this 
constraint is stronger than demanding that the entropy 
of the halo be larger than the entropy of the collection 
of particles that make up the galaxy halo in the early 
universe. This may be ascertained once the change in en- 
tropy is written as -k dfAD{f)/f IsqI where AD{f) 
is the change in the excess mass function. 



V. DARK MATTER HALOS AND THEIR 
CORES 

As an application of the above discussion consider a 
halo of DDM particles that has a King profile. These 
profiles have a core of almost constant density, a roughly 
l/r^ density run and finally a sharp drop in density to 
zero. The benefit of these profiles over the usual isother- 
mal profiles is that the King profiles have finite mass. 

The small galaxies in the local group have large Q val- 
ues. From a compilation of Local group dSph galaxies 
by Mateo infer that the largest value of Q is 

observed in Sculptor with a core radius of 110 pc and 
1-D stellar velocity dispersion of 6.6 km/s. Dalcanton 
and Hogan 33] show that one may assign a Q value of 
2 X 10"^ Mq/pc^/{ km s^'^ f to this galaxy. They also 
assign Q values to seven other local group galaxies with 
measured stellar velocity dispersions between 6 and 11 
km/s and My > —14. The Q values range between 10~^ 
and 10-''Mq/pcV( km s'^f. 

We choose a King profile for the dark matter halo with 
a concentration equal to 24 and central density set by the 
observed stellar profile core and stellar velocity dispersion 
[3^ . The total mass of the halo thus chosen is 2 x 10^ M© . 
We plot the excess mass function, normalized to the to- 
tal mass, for this profile in Figure ^ We also plot the 
primordial excess mass function that is greater than, but 
barely so, than the halo excess mass function. 

The King density profile drops to zero exponentially. 
The drop is sharp enough that for some small value of /, 
the King excess mass function will always be larger than 
the primordial one. Therefore, in the above comparison 
we only look at / values such that M(> /) < 0.99M(> 
0). Any realistic profile will have a more benign behavior 
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at small /; indeed, the excess mass function constraint 
will impose such behavior. 

The above exercise indicates that the decay parameters 
must be such that Qprim > 2x 10~^ Mq/pc^( km s~^)~^. 
We note that if the halo of Sculptor is less massive than 
we have assumed, then Qprim will have to be larger. If 
the halo is more massive, then Qprim can be lower. In 
scenarios with a small scale cut-off in the power spec- 
trum, the second possibility is more likely because the 
cut-off leads to a paucity of small mass halos 35j . 

We also note that requiring the excess mass function of 
the galaxy to be smaller than the fine-grained excess mass 
function in the early universe is a minimal constraint. 
The final effect on the galaxy of a small primordial phase 
space density is likely more complicated because of halo 
mergers. 

An interesting exercise at this point is to com- 
pare the constraint above to that obtained by de- 
manding that the entropy of the collection of par- 
ticles that make up the galaxy can only increase. 
We define entropy for N indistinguishable parti- 
cles as S' = -kN J d^xd^pf{p,x)[ln{f{p,x)h^) - 1]. 
For dark matter from decays, this works out to 
kN\n[27:e^^''^'^m'^Q~^h~^], which is very close to the 
entropy for an ideal gas of particles with the same 
Q. A comparison with the King profile entropy for 
Sculptor galaxy yields the result that in order for 
the entropy to not decrease we must have Qprim > 
10~^ Mq/ pc"^/( km s""'^)'^, weaker than the constraint 
obtained from the excess mass function analysis. 

The study by Dalcanton and Hogan [s^l found that 
the lowest mass halo core densities could be interpreted 
as resulting from a finite primordial phase space den- 
sity. However, they found no compelling argument to 
attribute to the core density vs. velocity dispersion they 
found over a wide dynamic range to the effect of a small 
primordial phase space density. 

The effect of a finite Q goes beyond just smoothing 
the inner cores of small galaxies. Larger galaxies are 
a result of the merger of smaller galaxies and hence the 
cumulative effects on the halo density of larger dark mat- 
ter halos could be substantial. Zentner and Bullock [s^ 
used semi-analytic arguments and the merger tree for- 
malism to show that a small phase space density can 
lead to a lowered concentration for dark matter halos. 
For the fermionic warm dark matter candidates, they 
calculate c(M) for a given virial speed and Q. We find 
that c(M) cx Q'-^-^^ approximates their results well for 
the range of Q values of interest. The exact relationship 
between c(M) and Q is more complicated in the present 
model as compared to the warm dark model case. This is 
because the constraint on Q depends on the final phase 
space profile of the halo. A more detailed study is re- 
quired on this subject. 

Why is this important? The concentration of ACDM 
halos as obtained from fits to low surface brightness 
galaxies (LSBs) falls on the low side. This might be due 
to two reasons. One, LSBs are extremely strongly bi- 




phase space density, f 



FIG. 1: The solid curve shows the excess mass function for 
a dark matter halo with a mass of 2 x 10^ Mq and a King 
profile. The dashed curve shows the excess mass function for 
the primordial phase space distribution of dark matter from 
decays with Q = 2 x 10"'^ Af0/pc^( km s"^)"^. 

ased to forming in the lowest concentration halos. Two, 
CDM is not the correct description of dark matter. As 
yet, we have no compelling model advocating that the 
first reason is correct. Clearly, the model presented here 
has the right features to help explain the observed low 
concentrations . 

The small scale end of the CDM power spectrum has to 
contend with another issue. CDM predicts a lot of sub- 
halos in the halos of the kind that would host the Milky 
Way. However, we see almost a factor of ten smaller 
number of galaxies .Sfj. On the other hand CDM sub- 
structure seems to be required to explain the strong lens- 
ing anomalies [sll . Plausible astrophysical solutions [s^ 
to this problem certainly exist, as do more exotic ones. 
Kamionkowski and Liddle advocated reducing small 
scale power to explain the lack of small galaxies in the 
halo of the Milky Way. Zentner and Bullock [s^ looked 
at this issue in detail and worked out the substructure 
fraction for linear power spectra with a small scale cut- 
off. Their semi-analytic results applied to our model of 
dark matter from early decays, and the results of the 
WDM simulations by Colin et al. suggest that this 
kind of dark matter will help alleviate the above men- 
tioned discrepancy. 

Another problem that has frequently been discussed 
is the observational case against cuspy halos that CDM 
simulations predict 42, 43, 44, 4^. Low surface bright- 
ness and dwarf galaxies seem to be better fit with a 
density profile that flattens towards the center (e.g.. 
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|43,|4^l49|). Such evidence for a core in higher mass 
galaxies is lacking. However, it is also true that that the 
larger galaxies are dominated by baryons in the centers 
making it harder to detect a core if it were present. Dark 
matter from decay of WIMPs can introduce a core in the 
smallest mass halos ~ 10^ Mq if Qprim is small enough, 
Qprim ~ 1O~*M0/ pc'^/( km s~^)'^, and reduce the con- 
centration parameter of larger halos. A model with such 
low Q values would produce a cut-off in the power spec- 
trum given by Ac = 0.25 Mpc. Is this consistent with 
having the universe reionize fully at a rcdshift of six? 
We explore this in the next section. 



VI. REIONIZATION 

Reionization provides a stringent constraint on the 
small scale power spectrum. The universe is reionized 
to a redshift of 6 and so there should at least be enough 
small scale power to do that. The recent WMAP results 
hint that the optical depth to Thomson scattering 
might be large (~ 0.1) indicating early reionzation. If 
true, this would have dramatic implications for the pa- 
rameter space of dark matter from decays. 

Barkana, Haiman and Ostrikcr 50] considered the con- 
straints on warm dark matter candidate from cosmologi- 
cal reionization. Their detailed considerations led to the 
conclusion that WDM with masses larger than about a 
keV could reionize the universe at redshift six. 

Here we perform a simple analysis to understand the 
effect of the decay on reionization using the semi-analytic 
models of Haiman and Holder [HJ • The ionized fraction 
is written as 



^Hii(^:) = Ph{z) / dz' 

J OO 

+ (1-Fhii(z'))^^^^(^') 



■ ) 



dz 



1/hii(z',z), (12) 



and Chii = 10 for the DM from decays. A robust calcu- 
lation will require input from n-body simulations. 

The efficiency parameter is the product of the average 
(over IMF) number of ionizing photons over the lifetime 
of the source, the escape fraction of ionizing photons and 
the fraction of baryons converted to sources (stars) . We 
take the efficiency to be 100. Given standard IMFs for 
metal-rich stars, we expect 4000 ionizing photons per 
baryon. Thus an efficiency of 100 implies escape fraction 
times fraction of baryons in stars of 2.5% in keeping with 
what is observed in the local universe. 

With the above inputs, the ACDM model gives a opti- 
cal depth of 0.091. For Ac — 0.1 Mpc//i, we get an optical 
depth of 0.061, while for Ac = 0.3 Mpc//i we get 0.035. 
However, for Ac = 1 Mpc//i we get an optical depth of 
0.017, inconsistent with an ionized 2 < 6 universe. 
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where pb is the average baryon density and Vhii(^', z^tM 
is the volume of region ionized at redshift z' in a halo of 
mass M. The subscripts la and lb 51] refer to different 
virial temperature ranges. Type la halos have tempera- 
tures between lO'' K and 2 x 10^ K, while lb halos have 
temperatures higher than 2 x 10^ K. In principle there is 
also a contribution from halos with virial temperatures 
smaller than 10^ K. This aspect has gotten a lot of atten- 
tion recently because it could be very important for an 
early reionization epoch. We are neglecting this contri- 
bution here to make a better comparison between CDM 
and DM from decays. In models with suppressed small 
scale power, this contribution is small. 

The reionization model as written above has two free 
parameters: Chii, e- Chii is the clumping factor for the 
IGM predicted by CDM simulations to be of order 10 at 
high redshifts. In the decaying DM scenario, the clump- 
ing factor may be smaller, thus impeding recombinations. 
For our simple analysis, we will take Chii = 15 for CDM 



FIG. 2: The curves show the collapsed fractions of halos with 
virial temperatures larger than 10* K. Below this tempera- 
ture, the gas can only cool by molecular hydrogen. The solid 
curve corresponds to a standard ACDM model with no tilt. 
The dashed curve corresponds to the same model but with 
a linear power spectrum modified according to Equation 1101 
with Ac = 0.1 Mpc//i. The dotted curve has Ac = 1 Mpc//i. 

The calculations of this section show that in order 
to be able to reionize the universe by redshift six, we 
expect that Ac < 0.3 Mpc/ft., which implies Qprim ^ 
10~^Mq/ pc"^/( km s^^)'^. The constraint is similar to 
the one we found in the last section by considering the 
phase space constraint from local group galaxies. If the 
universe is reionized at higher redshifts than six, this con- 
straint will tighten considerably. 

How reliable are the above calculations? Along with 
the simplistic modeling of the reionization process, one of 
the weak points of the above calculation is the calculation 
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of the collapsed fraction. We calculate cr^ (M) using the 
modified linear power spectrum and use the extended 
Press-Schechter scheme of Sheth and Tormen to 
calculate the mass function of dark matter halos. The 
results of this calculation for the collapsed fraction of 
halos with T > 10^ K are shown in FigureEl It is unclear 
that such a calculation will yield accurate results on small 
scales where the suppression in power becomes important 
and where structure formation is not fully hierarchical. 



VII. MIXED MODELS 

In scenarios of the kind we have discussed above, we 
may expect some of the DM to be produced during re- 
heating after inflation. Let us take the case of the grav- 
itino DM. If the reheating temperature is around 10^° 
GeV then we expect comparable contribution to the to- 
tal number density of gravitinos today from the decay 
and the reheating phase (see Roszkowski and Austri jg 
for recent work on early universe constraints on these 
models). The gravitinos from the reheating phase would 
behave as Cold Dark Matter particles. 

The main impact of these "mixed models" is to soften 
the small scale effects. This was considered in detail 
by Profumo et al. Q for the charged particle decay 
scenario. We write the dark matter density today as 
Pdm = /pDEC + (1 - /)PCDM where / is the fraction of 
the dark matter that results from the decay. The z — 
transfer function for the total dark matter is shown in 
Figure 01 Note that the suppression on small scales is 
reduced as expected. The form of the transfer function 
cannot be trivially obtained from the transfer function 
for the / = and / = 1 cases. 



VIII. DISCUSSION 

The phenomenology of DDM models is rich and much 
work remains to be done. To make robust contact with 
data on small scales we need large scale structure simu- 
lations and semi- analytic models. A few detailed simula- 
tions in this re gard have been run for models with keV 
WDM particle l25l l41l . with linear power spectrum cut- 
off on small scales |5J| , and to understand the dependence 
of the power spectrum on substructure . 

We looked at some beneficial features of these models 
in the previous sections and also the constraints on them. 
What are the other constraints? Narayanan et al. [s^ 
show that in the WDM scenario, masses smaller than 
750 eV are disfavored by the Ly-a data. Recent work by 
Viel et al. j57| peg the lower limit on the WDM mass at 
550 eV at 2a. If we naively compare the variance cr(M) 
in the WDM and the present model, then this may be 
turned into a constraint on the phase space density pa- 
rameter (see EquationlHJ), Q > 10~^ Mq/ pc'^/ ( km s"^)"^ 
corresponding to WDM masses larger than 750 eV, and 




k (h/Mpc) 



FIG. 3; The curves show the power spectrum for different 
values of /, the fraction of dark matter today that came from 
decays. The solid curve shows the / = case (CDM). The 
dashed curve shows the / = 0.5 case while the dotted curve 
shows the / = 1 case. It is clear that the suppression on small 
scales is much reduced for the / = 0.5 case. For comparison, 
we also plot (see thin solid curve) the transfer function for a 
1 keV Warm Dark Matter model. 



Q > 5 X 10^^ A/q/ pc'^/( km s~^)^ corresponding to 
WDM masses larger than 550 eV. 

Thus, we may summarize the present constraints on 
models where the bulk of dark matter today results 
from decays with a lifetime of about a year or smaller 
as Q > lO"** M©/ pc^/( km s~^)'^. Future comparisons 
with data will require large simulations to understand 
structure formation in these models; it does not proceed 
in a bottom- up hierarchical manner psl l4ll Is^ l . 

Can we push down to 100 kpc or even 10 kpc cut-off 
scales? Strong lensing flux anomaly technique is a sen- 
sitive probe of the amount of substructure [Ijl • Robust 
limits from this technique will require detailed theoretical 
understanding of issues such as the anisotropic distribu- 
tion of substructure in dark matter halos [H^ . Weak lens- 
ing constraints from future surveys is another promising 
avenue to learn more about dark matter on these small 
scales Q- 



IX. CONCLUSIONS 

In this paper, we have explored the important cosmo- 
logical consequences of dark matter from early decays - 
the cut-off in the power spectrum on small scales, and the 
limit on the phase space density of dark matter in halos. 
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The phenomenology is rich and there are multiple ways to 
search for the effect of such dark matter on structure for- 
mation. We pointed out that these models can suppress 
small scale substructure, create constant density cores in 
small mass halos due to the phase space constraint and 
reduce the concentration of larger mass halos. The mod- 
els that give rise to these early decays in supersymmetric 
theories are natural and they inherit the correct cosmo- 
logical abundance from the WIMPs that decay into them. 
Future observations of structure on small scales may be 
able to distinguish between cold dark matter and the 
dark matter from early decays. 
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